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Abstract Let N{L) be the number of eigenvalues, in an interval of length L, of 
a matrix chosen at random from the Gaussian Orthogonal, Unitary or Symplectic 
ensembles of A/" by A/" matrices, in the limit J\f oo. We prove that [N{L) — 



{N{L))]/ ^y log L has a Gaussian distribution when L oo. This theorem, which 
requires control of all the higher moments of the distribution, elucidates numerical 
and exact results on chaotic quantum systems and on the statistics of zeros of the 
Riemann zeta function. 
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Ensembles of N dimensional Gaussian Random Matrices (GRM) with invari- 
ances under the Orthogonal, Unitary or Symplectic groups corresponding to the 
GOE, GUE and GSE were introduced by Wigner and developed by Porter, Dyson, 
Mehta and others [1,2]. Wigner 's inspired surmise that the statistics of eigenval- 
ues of these GRM can be used to model the statistical properties of the observed 
spectra of complex nuclei turned out to be exactly right. There is indeed good 
agreement between the observed high energy level spacings, pair correlations and 
variance or A-statistics and those calculated analytically from the GRM in the 
limit A/" — > oo. Moreover, the GRM have been found to be the very robust "renor- 
malization group fixed points" of a large class of RM [3] which play an important 
role in many areas of physics and mathematics [1-7]. 

In the present work we focus on the large L (long wavelength) behavior of the 
random variable N{L) giving the number of eigenvalues of a GRM, chosen from any 
of the Gaussian ensembles, in an interval (y, y+L): We always consider the jV — > oo 
limit when the distribution is translation invariant and use units in which the mean 
spacing is unity. It is well known that the variance of N{L) grows like logL as 
L ^ oo. We prove that all the moments of ^(L) = [N(L) — L] / y/log L approach 
for large L those of a Gaussian distribution which implies (weak) convergence of 
^ (L) to a Gaussian random variable. We shall discuss later the connection of our 
result with the statistics of energy levels of quantum systems with generic chaotic 
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classical Hamiltonians and of the zeros of the Riemann zeta function [4,5,7]. 

It is a remarkable fact that the distribution of energy levels of the G(0,U,S)E 
are given by the Gibbs canonical distribution of the positions of charged point 
particles on the line interacting via the (two dimensional) repulsive logarithmic 
Coulomb potential, v{r) = — logr, at reciprocal temperatures /?= 1,2,4 respec- 
tively [1,2]. The particles with positions Xi, i = l,...,jV, on the real line, are 
confined by a uniform negative background, which produces a harmonic potential. 
The total potential energy of the system is 



In the corresponding circular ensembles of Dyson the Xi lie on a circle of length 
jV and the energy is given by the second term on the r.h.s of (1) with distance 
measured in the plane. The canonical Gibbs measures corresponding to Gaussian 
and circular ensembles become equivalent in the thermodynamic limit, jV" — > oo, 
yielding the same /c-point, /c = 1, 2, 3, .. correlation functions, for all /3 > [1,2,8]. 

These infinite volume correlation functions are known explicitly for the "solv- 
able" cases /? = 1, 2, 4 corresponding to the GRM. Defining as usual nj{xi, Xj) 
as the joint density for j-tuples, the corresponding Ursell functions 
Uk{xi,...,Xk) [9] are given by 




(1) 



Ui{xi) = ni{xi) = 1 



(2) 
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U2{xi, X2) = n2{xi, X2) - ni{xi)ni{x2) = n2{xi - X2) - 1 



t/fc(a;i,...,a;fc) ^ J](-ir-i(m- l)!n^G,({a;G,}) 

G 3=1 

where G is a partition of the indices (1, 2, .., A;) into m subgroups Gi, .., Gm and 
{xGj} are the xi with indices in Gj. (In the GRM literature, Yfc = (— 
is usually called the k^^ cluster function). The integrals, Uk, of Uk over a k- 
dimensional cube having sides of length L are directly related to the cumulants 
Cj{L),j — 1, ...,k of the random variable N{L), the number of points (or eigen- 
values) in an interval of length L, which we shall take for definiteness to be the 
interval [-t,t], L = 2t. Thus 

Ui = J Ui{xi)dxi = 2t= {N{L)) = Ci{L) (3) 

U2= [ [ dxidx2U2{xi,X2) = {N{N-l))-{Nf = {{N-{N)f)-{N) 
J -t J -t 

= C2(L)-Ci(L) 
Us = CsiL) - 3C2(L) + 2Ci(L) , etc 

Using the generating function [9] 

00 

F{^^) = Y,En{L)e-', (4) 

n=0 

with En{L) the probability of having exactly n particles in the interval L, we have, 

00 „ 00 ^ 

n=l n=l 
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This gives 



Ck = J2^k,jCj + {k - l)!Ci + tJk 




(6) 



where 



b. 



7 = b, 



{n-l)b, 



2<j<n-l 



b. 



n,n 



= -1 



n > 2 



The Uk take on a particularly simple form for the GUE (/? = 2) [1,2], 

k 

Perm^=i 

where S{x) = (tto;)"^ sinTra; and Xk+i = xi, so the indices are to be thought of as 
being on a circle. Using (6) one readily obtains [1,2] 



Note that for a system with short range interactions C2{L) would grow like L 
but the logarithmic interactions between the (charged) particles induce a sort of 
local crystalline order reducing the variance to logL. It is this strong correlation 
which produces level repulsion between the eigenvalues and makes the large scale 
behavior of the fluctuations far from obvious. Defining now the normalized random 
variable 



C2(L) = (logL)/7r2 + 0(l). 



(8) 
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the k^^ cumulant of r]{L) will be = Cfe/[logL/7r^] 2 . As is well known c.f. [10], 
r]{L) will approach a Gaussian random variable (with mean zero and unit variance) 
as L — > 00, if and only if all Ck, k > 3, go to zero, i.e. if Ck = o([logL/7r ja). 
Using an induction argument based on the recurrence relation (6) and the equality 
(7), this corresponds to proving that 

Sk{t) = / ... / dxi..dxkS{x2 - xi)S{xz - X2)...S{xi - Xk) 

J-t J-t (10) 

= 2t + o[{logt)^, k > 3 
We shall actually prove that Sk{t) = 2t + 0{\ogt) which implies that, for A; > 3 
Ck{L) = O(logL); in fact we believe that for A; > 3, Ck{L) stays bounded as 
L — i> 00, as suggested by the explicit asymptotic evaluation of the integrals. 

53 = 2t--^logt + 0(l) 

(11) 

54 = 2t- ^\ogt + 0{l) 
which gives, using (6), that C3 and C4 are of 0(1). 

To prove (10) we make use of the fact [1] that 

Sk{t) = J2>^iit)=T^^'(t) (12) 

where the Xi{t) are the eigenvalues of the integral operator 
S(t), (S/)(a;) = J^^dyS{x — y)f{y). It is known [1] (and can be easily proven) 
that the spectrum of S lies in [0, 1]. We can now use an induction argument to 
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prove that Sk{t) = 2t + 0{logt). This is so for A; = 1, 2 (also for 3, 4), and for A; > 2 
we have 

TrS^+^ = TrS''-Tr[S^-^(S-S2)]. (13) 

The first term is of the desired form by the induction assumption while the terms 
in the parenthesis are positive operators, and ||S^~^|| < 1 so it can be taken out 
of the product yielding Tr S'^"''-'^ = Tr S*^ + O(logt) and the proof is complete. 

Our results readily extend to show that if we divide up the real line into a union 
of intervals of length L, let Nj{L) be the particle number in [jL, (j + 1)L] and set 
r]j(L) = [Nj(L)—L]/y^C2{L), j ^ Z then the {r]j(L)} approach, as L — > oo, jointly 
Gaussian random variables with mean zero and covariances {rjjTjk) = Sj^k — |<^i±i,fc- 

To prove the results for the GOE, /3 = 1, we use an identity conjectured by 
Dyson and proved by Gunson [11] (we thank Freeman Dyson for pointing this 
out to us) . According to this identity, superimposing two noninteracting Coulomb 
gases, say blue and red, in the circular ensemble at reciprocal temperature P = 1 
and then looking only at alternate particles, e.g. at all the odd (or even) ones, 
yields the distribution at /3 = 2. Considering now the number of particles in an 
interval of length L gives 

NtotaiiL) = Nlll{L)+Ni.l\{L) = N^I\{L) + N(IUl) = 2N^'\L) + {0, ±1) (14) 
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where the superscripts (1,2) stand for the random variables obtained from the 
ensembles with (3 = 1,2, and Ni,iuf,, A^r-ed ^ire independent. This shows immediately 
that in the infinite jV limit, the variables N^'^\L) and N^^\L) normalized by 
the square root of their variances have the same asymptotic behavior. Taking 
{N^Lm = {n!:11{L)) = (Ar(i)(L)) = L we have 

kiKliL) -L) + (ivW (L) - L)f) = ((iV(i)(L) - Lf) = 

(15) 

2((Ar(2)(L)-L)2)~^^ 

giving the well known variance of the GOE [1] . 

For the GSE, /3 = 4, we use the equality between the statistics of its eigen- 
values and the odd eigenvalues of the GOE [1]. This again leads to Gaussian 
asymptotics with a variance given by log L + 0(1). It seems very reasonable to 
expect and one can give strong heuristic arguments, based on the "long wavelength 
response" of "Coulomb" systems, that the Gaussian nature of the fluctuations, 
with variances (2/7r^/9) log L, holds for all (3. ( We thank Bernard Jancovici for 
pointing this out to us, see also [8].) 

Using more detailed information on the spectrum of S (see [12]), it follows that 
Sk{t) = 2t— 7r~^ Sj=i logt+o(logt) which using (6) implies that Ck = o(logt) 
for all A; > 3. (We are indebted to Harold Widom for this information). Widom 
also noted that our proof does not make any use of the specific form of S. It only 
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uses the property spec(5') e [0, 1] and the fact that Tr(S — S^) —>■ oo (as t — > oo). 
The conclusion therefore holds for a larger class of matrix models [3] . 

The local statistics of the eigenvalues ej,j = 1, 2, ...,oo, < ei < €2 < ... of 
a classically chaotic quantum Hamiltonian (CQH), such as the geodesic flow on a 
(non-arithmetic) surface of constant negative curvature or the Sinai billiard appear 
to coincide at high energies with those obtained from the GRM [4] . More precisely, 
if we consider the energy levels of a generic CQH, suitably scaled so that the mean 
distance between levels is unity, in an interval {y, y + L) then their statistics, 
obtained by letting y vary uniformly in some interval (Ti,T2) will coincide, for 
T2 ^ 00, with that obtained from one of the standard GRM ensembles when the 
matrix size tends to infinity. Our result then predicts a Gaussian distribution of 
the fluctuations in the number of levels N[L) when L — > 00. A numerical check of 
this for some CQH will require using energy levels in a (scaled) energy interval L 
with 1 « L « (T2 - Ti) , Ti large. 

It is also interesting , as emphasized by Berry [4], to consider in addition to 
the local statistics of quantum levels also their global statistics. These correspond, 
in our context, to fluctuations in the number of levels in an interval {y, x) whose 
length L{x) — x — y{x) is not fixed but grows with a; as a; varies in the interval 
(Ti,T2) with T2 — > 00. This includes in particular the case y = 0, L{x) = x 
corresponding to the fluctuations in the number of eigenvalues less than x. This 
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quantity, normalized by the square root of its variance, was conjectured in [5] 
(where it is denoted by Nfi{x)) to have a Gaussian distribution as a; ^ oo for all 
CQH. If true this would be a general characterization of CQH and distinguish them 
from integrable systems where it was found rigorously that the global distribution 
is non-Gaussian [13]. Quite generally, it was shown by Berry [5] that when L{x) > 
-^max(a;), the variance of N{L{x)) saturates for L > Lmax- Berry also found that 
Nji{x) (averaged over some interval containing many eigenvalues but very small 
compared to x) grows for billiard systems like (27r)~^ logo;. 

As already noted, the distribution of eigenvalues in the GRM is translation 
invariant, when the matrix size jV goes to infinity, so there is no analog of Lmax in 
our considerations. One can however consider fluctuations in N{L) for an interval 
L{J\f) which contains a number of eigenvalues small compared to J\f but goes to in- 
finity when jV" — > oo, e.g. in the circular ensemble we could have L{N) ~ N'^, 7 < 1 
or even like aj\f, a << 1. For the Coulomb system with neutralizing background 
it is also possible to consider semi-infinite systems with various boundary condi- 
tions and/or nonuniform background. Some such systems have been considered in 
[8] and we believe that our results about Gaussian behavior would extend also to 
these systems which might model some of the saturation features of CQH. 

We turn finally to the (non trivial) zeros of the Riemann zeta function ((z) = 
'^^^i Ti~^ which are, according to the Riemann Hypothesis, of the form = \ + i'ji 
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As pointed out by Berry [4] there are reasons to expect similarities between the 
statistics of the 7„ and of energy levels of CQH. In fact Montgomery [14] proved 
that the pair correlation function of the 7^ agrees with that of the GUE, Eq. 
(7). Numerical calculations by Odlyzko [7] give striking evidence that the nearest 
neighbor level spacing distribution of the 7^ is, for large n, indeed the same as 
that obtained from the GUE. In a very interesting recent paper Rudnick and Sar- 
nak [7] greatly extended the results of Montgomery by showing that the n-point 
correlation functions of these zeros converge, on a large class of test functions to 
those of the GUE. Moreover the normalized global fluctuation in these zeros cor- 
responding to Nfi{x), was shown by Selberg [15] to have a Gaussian distribution. 
The same arguments imply that the local fluctuation in their number in an inter- 
val {y,y + L) averaged over (Ti,T2) and properly scaled becomes Gaussian when 
T2 — > 00, followed by L — > 00 (we are indebted to Peter Sarnak for explaining 
this to us). The results proved here thus fit completely with the picture of the 
statistics of Riemann zeros being in the same "universality class" as that of the 
GUE. 
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